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FOREWORD 


This report constitutes the final part of the work completed on the 
research project "Determination of Atmospheric ?( llutants from Infrared 
Radiation Measurements." The work was supported by the NASA/Langley 
Research Center through research grant NSG 1282. The grant was 
monitored by Mr. John T. Suttles of the Experiment Analysis Branch 
of the Atmospheric Environmental Sciences Division. 
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GROUND TEMPERATURE MEASUREMENT BY 


PRT-5 FOR MAPS EXPERIMENT 
By 

S.K. Gupta^ and S.N. Tlwarl^ 


SUMMARY 

A simple algorithm and computer pr^'-gram have been developed for deter- 
mining the actual surface temperature from the effective brightness temperature 
as measured remotely by a radiation thermometer called IRT-S. This procedure 
allows the computation of atmospheric correction to the effective brightness 
temperature without performing detailed radiative transfer calculations. 

Model radiative transfer calculations are performed to compute atmospheric 
corrections fo: several values of the surface and atmospheric parasnters 
Individually and In combination. Polynomial regressions were performed 
between the magnitudes or deviations of these parameters and the corresponding 
computed corrections to establish simple analytical relations between them. 
Analytical relations have also been developed to represent combined correction 
for simultaneous variation of parameters In terms of their Individual 
corrections. These analytical relations can be used In all future work 
to determine the actual surface temperature without having to perform 
expensive radlati^ve transfer calculations. Model calculations have also 
been made to examine the sensitivity of the retrieved surface temperatures 
to uncertainties in the various surface and atmospheric parameters. 

^ Research Assoc..ate, Old Dominion University Research Foundation, Norfolk, 
Virginia 23508. 

- Professor, Department of Mechanical Engineering and Mechanics, Old Dominion 
University, Norfolk, Virginia 23508. 


1. INTRODCCTI(»l 


McMurwBcnt of conctntrationa of ga««ou« «tTO«ph«ric pollutant! Ilka 
CO, CH^i MBj and SO 2 la of graat Inportanca bacauaa of thair conaldarabla 
affacr. on varioua natural procaaaaa affactlng human haalth, anvlronmant 
and aataorology [1]*. Succaaaful attampta hava baan mada racantly, to 
naaaura aavaral of thaaa pollutanta by paaaiva ramota-aanalng Inatrumanta 
nountad on aircraft [2,3] with tha Intantlon of avantually adapting than to 
aatalllta uia. Thaaa InatnoMnta maaaura tha upwalllng Infrarad radiation 
in an appropriata apactral ragion, which in turn carriaa tha aignatura of tha 
particular pollutant aolaculaa [4]. 

Extraction of tha pollutant concantrationa from thaaa maaaurad aignals 
i.'aqulraa, in addition to aophlaticatad data raduction algorlthma [5], prior 
luiowladga of aavaral paraiMtara of tha atmoaphara aa wall aa tha undarlying 
aurfaca. Tha atooapharic paramatars, a.g., praaaura, temparatura /*.na humidity 
profilaa can ba obtalnad althar from additional inatrumanta mountad aboard 
tha aircraft Itaelf or from othar indepandant aourcaa (a.g., radloaondaa) . 

Tha aurfaca paramatara, namaly, tha emittanca and tamparatura ara not avall- 
abla from any auch aourcaa. Any maaningful phyaical contact maaauramant of 
Chesa paramatara la not faaaibla bacauaa of the vaat areaa covered in remote 
aanalng exparimanta and large apatlal variatlona of tha quantitiaa Involved. 
Reliable eatlmatea of aurface emittanca, however, can be obtained from an 
approximate knowledge of the color, texture and compoaitlon of the aurface. 
Thla ia poaalble becauae extenalve tabulation relating baalc characteriatlca 
of the aurfacea to their emittanca values are available in the literature [6]. 
Thla leavea the aurfaca temperature aa the important unknown quantity to ba 
determined . 

'^Nundjera in bracketa indicate referencea. 


■) 


The work pr«s*nt«d in thli report wet ptrforaad in support of the 
project MAPS (HeeeurwBMint of Air Pollution from Setellitee) et the Lengley 
Reaeerch Center. Current work under this project involves frequent flights 
of e ges-filter correlation radiometer (CFCR) on an aircraft. Thia instru- 
ment is designed to measure carbon monoxide concentration in the troposphere. 

A radiation thermometer called PRT-S (Precision Radiation Thermometer-5, 

Barnes Engineering Company, Stamford, Connecticut) is flown alongside the 
GFCR for measuring the temperature of the underlying surface. This latter 
instrument is essentially a broad-band radiometer which operates in the 
11 micron atmospheric window. Under certain ideal conditions, namely (i) when 
the emittance of the target surface is unity, and (11) when the attenuation 
by the atmosphere between the target surface and the sensor is negligible, 
this instrument .an direccly read the surface temperature. 

These ideal conditions, however, are always far from being met in the 
real world situations. In addition to the deviation of surface emittan:e 
from unity, line and continuum absorption by water vapor in this region [7,8] 
introduce error in the temperature measurement by PRT-5, which needs to be 
determined and corrected for. A strong temperature dependence of the continuum 
absorption coefficient [8,9] introduces additional uncertainty in this measure- 
ment. Andlng and Rauth [10] and Prabhakara et al. [11] have discussed the 
use of differential absorption by water vapor in adjacent spectral regions 
to correct for the attenuation due to water vapor. Since PRT-5 is a single- 
channel, non-dlsperslve Instrument, the above technique is not applicable in 
the present case. 

The purpose of the present study is to develop algorithms and software 
for computing correction to the measured surface temperature arising from 
the above-mentioned causes. The attenuation due to water vapor absorption 
may be evaluated accurately using a line-by-line radiative transfer model. 
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Routina data raduction uaing llna-by-lina nodal canda to ba axtranoly 
cxpanalva [12] and, tharafora, inaxpanalva algor! thna hava to ba davalopad 
to dataxnlna tha nagnltuda of tha corraction. Thia can ba dona by parfomlng 
aona nodal c.. ' ulationa and uaing thoaa raaulta to aatabllah alapla analytical 
ralationa batwaan tha nagnikudaa and daviationa of tha paranatara fron thair 
atandard valuaa and tha corraaponding corractiona raquirad to coaipanaata for 
than. Thaaa ralationa can than ba uaad in tha futura to conputa tha corrac- 
tiona for uiy aat of paranatara without parfoming datailad radiativa tranafar 
calculations. 

Ralavant thaory and tha davalopnant of algorithm is prasantad in 
Saction 2. Data sourcas and computation procedure are referred to in 
Section 3. Results of tha nodal calculations are discussed in Sections 4 and 
5. Application of the algorithm is presented in Section 6. Section 7 
contains results of some sensitivity calculations. Important conclusions 
of this work are discussed in Section 8. 


2. THEORY AND ALGORITHM 


Th« qu«nticy maturtd by any broad-band radlosMtar la tha Intagractd, 
upwalling radlanca which for a hovoganaoua nonacattaring atmoaphara can be 
axpraaaad aa [4,13] 

U)2 

E ■ I (2.1) 

whara uii and U 2 rapraaant tha fraquancy llmlta of tha obaarvad apactral 
ragloQ and f(u) la tha inatruaant filtar function. Eq(u) tapraaanta tha 
aonochromatlc thamal rr.dlanca amlttad by tha underlying aurfaca and tha 
atmoaphara and la given by 

Eg(w) - c(w) B(o)»T^) t(u),0) 
h 

+ j B[w,T(z)) fdT(u),z)/dz] dz (2.2) 

*^o 

where C(u)) la the surface emlttancc, B(u>,T) Is the Planck'a blackbody 
radlanca for temperature T , T^ la the aurfaca temperature, T(z) la 
tha temperature at altitude z, and t(u,z) la tha monochromatic trans- 
mittance of tha atmosphere. The first term on the right-hand aide of 
equation (2.2) represents the radiation from the surface while the second 
term represents the radiation originating In the atmosphere. Solar radia- 
tion reflected from the surface and the scattered radiation make negligible 
contribution to the upwelling radiance In the 11 micron region. 

The quantity measured by a direct-reading instrument like the PRT-S 
is the effective brightness temperature (EBT) of the surface, denoted by 

T and la defined as 

e 

B(uj,T^) f(w) du) - E (2.3) 
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where the blade-body radiance corresponding to equals the total upwell- 
Ing radiance. Simple Inversion of equation (2.3) for radlometrlcally 
measured values of E would yield for any given altitude and set of 

surface and atm<' :'pherlc conditions. Inversion of equation (2.1) to extract 
the actual surface temperature on the other hand, would Involve 

detailed radiative transfer calculations. Furthermore, the accuracy require- 
ments of the atmospheric problems necessitate the use of llne-by-llne model 
for transmittance computation [lA]. 

2.1 Analytical Formulation 

The need to perform llne-by-llne radiative transfer calculations for 
every set of surface and atmospheric conditions may be eliminated by per- 
forming some model calculations and establishing relations between T^, 
and the magnitudes or deviations of the Input (surface as well as atmospheric) 
parameters. Since Is equal to for surface emlttance equal to 

unity and no water vapor In the atmosphere, these are adopted as the Initial 
conditions. The quantities which give rise to the difference between T^ 

and T (l.e., the correction denotec' by AT) are the magnitude of water 
s 

vapor burden and the deviation of surface emlttance (e) from unity. 

Analytical relations will be established, therefore, between the above 
quantities and the resulting corrections AT . The correction may be expressed 
as 

AT - T - T (2.4) 

e s 

where T^ is obtained by inverting equation (2.3). Since radlometrically 
measured values of E cannot be obtained for all the desired conditions 
for establishing the necessary relations, they are computed from equation (2.1) 
for a series of values of the surface temperature, T^ and several values of 
each of the input parameters. 
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Gupta and Tiwarl [15] analyzed the magnitudes of corrections to EBT 
caused by deviations of various input parameters for radiometric measurements 
made In the 4.66 micron region and found that the correction can 
usually be expressed In terms of the deviation Ax (of the parameter x) 
by simple analytical functions. For most cases, these relationships may 
be represented by a polynomial as 

N 

AT - 23 a^(Ax)" (2.5) 

n"o 

where N Is a positive integer. The first few terms of this polynomial 
should adequately represent the correction AT for the expected range of 
the deviations Ax . The coefficients can be obtained by polynomial 

regression in the results obtained from model calculations. Future compu- 
tations of AT can then be made from a^ and Ax without having to perform 
the lengthy radiative transfer calculations. Results of the model calcula- 
tions for the deviation of a single parameter and their analysis Is presented 
In Sec. 4. 

A more complicated situation arises when more than one Input parameters 
undergo simultaneous deviations. The combined correction, to be denoted by 
AT^ may be represented as a function of the individual corrections. For 
example, if Ax and Ay represent deviations of two input parameters x and 
y from their initial values, the combined correction may be expressed as 

AT, - f(AT , AT ) (2.6) 

0 X y 

where AT and AT represent the individual corrections caused bv the 
deviations Ax and Ay separately. Additional model calculations have 
been performed to examine the nature of the above function for different 
sets of parameters. Results of these calculations and their analysis are 
presented in Sec. 5. 


2.2 Water Vapor Absorption 

The total absorption coefficient k^(u) for water vapor in the 
11 micron window may be expressed as [9,11] 

k^(o)) - kj^(o)) + kp((D) (p-Ph^q) + k^(a (2.7) 


where 


k^Cui) - absorption coefficient due to water vapor lines, 
k^(u>) - continuum absorption coefficient due to pressure broadening 
by other gases. 


k^(u) - continuum absorption coefficient due to self-broadening, and 

p,p„ - the total pressure and partial pressure of water vapor respectively. 
n20 


Absorptions represented by k^ and k^ are also called the p-type and 
e-type absorptions respectively. Measurements of k^ and k^ are reported 
by several workers [7,8]. Roberts et al. [9] reviewed these and other recent 
laboratory measurements by Burch and reanalyzed a bulk of field measurements 
[16] to conclude that in the vicinity of 300®K, k^ is very small compared 
to k 


The total absorption may, therefore, be expressed as 

k^(u)) - kj^(o)) + k^(0)) Pjj^Q 


( 2 . 8 ) 


A strong negative temperature dependence of k^ was observed by Blgnell [8] 
and Burch [7]. Roberts et al. [9] have analyzed the above temperature 
dependence data and that from several <}ther workers and have shoim that 

kg(o)) - k°(o)) exp[T^(i - 2 I 6 )] (2.9) 

gives the value of k^(o)) at the temperature T , where k°(aj) refers to 
a temperature of 296 ®K and the best estimate of T^ is found to be 1800®K. 

The frequency variation of k°(ui) in the 8-12 micron region was also 
studied by Roberts et al. from an analysis of the experimental data from 
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the various available sources, Roberts et al. [9] have developed an 
empirical represent;.tion which can be expressed as 

kg(w) ■ a + b exp(-Soj) (2.10) 

where the constants a, b and 8 are 

a ■ 1.25 X 10“*^ mol“^ cm* atm“*, 
b ■ 2.34 X 10“^* mol“' cm* atm'*, 

and 

■ 8.30 X 10“* cm. 

It has been pointed out by the authors that while equation (2.10) can 
be used for the entire region up to 30 microns, the temperature dependence 
expressed by equation (2.9) may not be valid far beyond the 8-12 micron 
region. 


3. DATA SOURCES AND COMPUTATION PROCEDURE 


The model atmosphere adopted for the present work Is essentially the 
U.S. Standard Atmosphere, 1962 and the water vapor profile is approximately 
an average of winter and summer mid-latitude distributions [17]. Since 
other minor atmospheric constituents (e.g., CO 21 N 2 O) are inactive in this 
spectral region, they were excluded from consideration in the model atmosphere 
as well as the radiative transfer. 

For the present work performed to support aircraft observations, tropo- 
sphere up to an altitude of 17,500 ft. was considered and was divided into 
ten layers of unequal thicknesses. Pressure, temperature and water vapor 
concentration for these layers are given in Table 3.1. The altitudes at the 
layer boundaries are listed in the first column of the table. Data in the 
next three columns represent the averages for layers between the altitudes 
Just below and above these values. 

Table 3.1 


Model Atmosphere used for the Study. 


Altitude 

(ft) 

Pressure 
(mm Hg) 

Temperature 

(“K) 

Water Vapor Cone. 
(10^ ppmV) 

0 

763.16 

287.66 

7.626 

500 

742.65 

286.17 

7.240 

1,500 

716.03 

284.19 

6.726 

2,500 

677.55 

281.22 

5.972 

4,500 

628.87 

277.26 

5.100 

6,500 

583.07 

273.30 

4.222 

8,500 

540.03 

269.34 

3.338 

10,500 

499.59 

265.38 

2.669 

12,500 

461.66 

261.42 

2.074 

14,500 

417.57 

256.48 

1.493 

17,500 
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Spectral parameters of individual lines of water vapor in the 9.8- 
13.2 micron (1020-760 cm"’) region were obtained from the AFCRL line 
parameter compilation [18], and were used in evaluating the line component 
of the absorption coefficient using line-by-line model. The information 
presented in Sec. 2 was utilized to evaluate the continuum component 
(self-broadened) of the absorption coefficient. 

Integrated upwelling radiance was computed using equation (2.1). 

Radiative transfer calculations were performed using total absorption 
coefficient obtained by combining the line and continuum conq>onents. 

Effective brightness temperature is obtained by inversion of equation (2.3) 
using quadratic interpolation. The Instrument filter function as supplied 
by the manufacturer is shown in Fig. 3.1 and was used in radiance calcula- 
tion after digitization. A compu':er program named PRTFIVE has been 
developed to accomplish the above objectives and is listed in Appendix B. 

It evaluates in a single run, the Integrated upwelling radiance and therefrom 
the EBT for eight values of actual surface temperature from 290 to 325°K and 
for ten different altitudes at the top of each of the ten layers. Program 
runs were made for several values of each of the parameters covering realistic 
ranges for their variations. Analytical relations were developed by carrying 
out regression in the tables of Ax and corresponding tables of AT for 
the various parameters. Starting from a linear fit, higher degree fits were 
attempted if the magnitudes of the residuals between the computed values and 
those obtained from the fit exceeded O.l^K. 
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PERCENTAGE 




4. MODEL CALCULATIONS FOR ONE PARAMETER 


The Input parameters considered Important for this study are (1) 
surface emlttance, (11) water vapor distribution, and (111) atmospheric 
temperature profile. Model calculations have been performed and analyzed 
for deviation of a single parameter at a time and also for simultaneous 
deviations of more than one parameter. Results obtained for deviations of 
a single parameter are discusser xxl this section. Discussion of the results 
obtained for simultaneous variation of two parameters is covered in the 
next section. The corrections for each value of the parameters are obtained 
for eight values of the actual surface temperature, T^ and for all ten 
altitudes. 

4.1 Surface Emlttance 

Integrated upwelllng radiance and EBT were evaluated for five values of 
surface emlttance between 1.00 and 0.80. Computations have been made In 
absence of any water vapor In the atowsphere as well as In presence of a 
standard water vapor burden. The dry atmosphere results are presented first 
and Table 4.1 lists the values of EBT and corrections for all eight values 
of Tg for this case. Since there Is no attenuation of radiation in the 
atmosphere In this case, these values of EBT and correction are Independent 
of the altitude. Figure 4.1 shows the magnitudes of the correction as a 
function of the deviation of emlttance for T^ <■ 300, 310, and IZO^K. 
Polynomial regression In these data and those for other values of T^ shows 
that these results can be adequately represented by quadratic expressions. 

The coefficients of quadratic regression, a^, and a^ were determined for 
all values of T^ . It can be seen from figure 4.1 that all these curves 
pass through the origin and hence, the coefficient a^ Is Insignificantly 
small In all cases. The significant coefficients a^ and a^ for all values 
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of art prtttnttd In Tablt 4.2 and It ahowt that both cotfficlanta art 
dtptndtnt on tht actual turf act ttmptraturt, T . Figure 4.2 dtpicta tht 

B 

variation of a^ and with actual surface ttaptraturt. 


Tablt 4.2 



Bagrestion cotfflcitnta a^ and a 2 at functions of 
the actual surface temperature for the emittance vari- 
ation study in a dry atmosphere. 


Actual Surface 
Temp. (*K) 



Computations in the presence of a standard water vapor burden (or wet 
atmosphere) were also made for all eight values of T^ and all ten altitudes. 
Because of attenuation due to atmospheric water vapor, the EBT changes with 
the altitude in this case. Values of EBT and the correction for all values 

of T and the altitude of 17,500 ft. are listed in Table 4.3. Oiffeiences 

s 

between T^ and the values of T^ which correspond to c ■ 1.00 are due to 
the attenuation of the surface radiation by the atmospheric water vapor. 
Corrections due to deviation of emittance are, therefore, obtained by sub- 
tracting the values In the second column from those in columns 3-6 of this 


table and are listed right under the T^ values. 
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Tab It 4.3 


Efftctivt brighcntat camptraturt and corrtctlon ac 
a function of aurfact tmittanct for an altitude of 
17,500 ft. 


Actual 

Surface 

Temp.i’K) 

Effective Brlghtnoee Teaperatux 

• CK) 


e - 1.00 

c - 0.95 

e ■ 0.90 

e - 0.85 

c ■ 0.80 

290 

288.501 

285.600 

282.617 

279.556* 

276.415* 


0.000 

-2.901 

-5.884 

'8.945 

-12.086 

295 

292.910 

289.898 

286.805 

283.620 

280.343 


0.000 

-3.012 

-6.105 

-9.290 

-12.567 

300 

297.341 

294.214 

291.006 

287.707 

284. 306 


0.000 

-3.127 

-6.335 

-9.634 

-13.035 

305 

301.786 

298.546 

295.223 

291.806 

288.282 


0.000 

-3.240 

-6.563 

-9.980 

-13.504 

310 

306.246 

302.893 

299.453 

295.917 

292.275 


0.000 

-3.353 

-6.793 

-10.329 

-13.971 

315 

310.719 

307.255 

303.696 

300.040 

296.275 


0.000 

-3.464 

-7.023 

-10.679 

-14.444 

320 

315.205 

311.625 

307.949 

304.172 

300.283 


0.000 

-3.580 

-7.256 

-11.033 

-14.922 

325 

319.700 

316.006 

312.214 

308.313 

304.300 


0.000 

-3.694 

-7.486 

-11.387 

-15.400 




*Value« obtained by extrapolation 












Figure 4.3 thowr the veluee of the cempereture correction ee e function 
of the deviation of aurfece emittance for the altitude of 17|500 ft. end 
three veluee of 300, 310 end 320*K. Polynoaiel regreeeion in theee 
date ehrve again the curves can be represented by quadratic exprescions with 
residues ouch snaller than the acceptablt: 0.1*K. The regression coefficients 
a , a, and a> were determined for all values of T. and the altitxida and 

O 1 4 9 

all e^ are insignificant for the reason mentioned earlier. The coefficients 
a^ and a 2 are preaeated in Table 4.4. This table shows again that both 
coefficients a^ and a 2 ere dependent on the surface temperature, T^ as 
well as the altitude. As illustrations, the variation of these coefficients 
with T^ for the altitude of 17,500 ft. is shown in Figure 4.4. Figure 4.5 
shows the variation of these coefficients with altitude for T^ ■ 300*K. 

It is Important to note chat the emittance correction determined here 
for Che wet atmosphere case has to be combined with the correction due to 
the presence of water vapor. As a result, the total correction is this case 
is given by 

AT - ATg + AT^^ (4.1) 

where AT^ is the correction determined from the coefficients given in 
Table 4.4 and A7^ corresponds to the presence of standard water vapor 
burden. 

4.2 Water Vapor Distribution 

The water vapor is an important absorbing coustlCuent for infrared 
radiation in the troposphere and is the only absorber in the spectral region 
of this Instrument. Its variation, therefore, has significant effect on 
Che correction to EBT and is considered in two ways. These are (1, Increasing 
or decreasing water vapor concentration as multiples of the standard water 
vapor profile in all layers, and (11) increasing concentration in some layers 


19 



DEVIATION OF EMITTANCE 


Fig. 4.3 Variation of correction with the deviation of emittance for 
the wet atmosphere from 17,500 ft. 
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ACTUAL SURFACE TEMPERATURE, TsOO 



Fig. A. 5 Variation of regression coefficients a^ and with 
altitude for the wet atmosphere and Ts * 300*K. 


Table 4 .4 


Regression coefficients 


and a^ are functions of actual 


surface temi)erature and altitude for surface emlttance study. 


Altitude 


Regression Coefficients 


vtt.; 

290*K 

295*K 

300 “K 

305 "K 

310 *K 

315-K 

320*K 

325 “K 

500 

64.056 

66.117 

68.352 

70.563 

72.841 

75.147 

77.461 

79.859 


-20.357 

-21.916 

-22.333 

-23.044 

-23.551 

-24.047 

-24.610 

-24.871 

1500 

62.404 

64.484 

66.685 

68.935 

71.209 

73.487 

75.836 

78.216 


-19.302 

-20.802 

-21.441 

-21.996 

-22.536 

-23.144 

-23.476 

-23.872 

2500 

61.196 

63.286 

65.488 

67.735 

69.981 

72.277 

74.620 

76.994 


-18.541 

-20.011 

-20.671 

-21.248 

-21.897 

-22.429 

-22.769 

-23.164 

4500 

59.572 

61.625 

63.865 

66.107 

68.345 

70.694 

72.997 

75.316 


-17.525 

-19.189 

-19.670 

-20.269 

-20.934 

-21.146 

-21.695 

-22.252 

6500 

58.568 

60.619 

62.853 

65.062 

67.362 

69.660 

71.949 

74.244 


-16.894 

-18.571 

-19.074 

-19.778 

-20.191 

-20.564 

-21.139 

-21.787 

8500 

57.961 

60.007 

62.235 

64.485 

66.734 

69.024 

71.292 

73.586 


-16.505 

-18.203 

-18.725 

-19.262 

-19.843 

-20.227 

-20.891 

-21.480 

10500 

57.617 

59.660 

61.931 

64.129 

66.399 

68.647 

70.918 

73.176 


-16.276 

-17.995 

-18.353 

-19.070 

-19.457 

-20.122 

-20.718 

-21.458 

12500 

57.418 

59.458 

61.725 

63.920 

66.186 

68.430 

70.697 

72.949 


-16.137 

-17.874 

-18.244 

-18.963 

-19.349 

-20.022 

-20.623 

-21.367 

14500 

57.309 

59.346 

61.611 

63.803 

66.067 

68.308 

70.572 

72.821 


-16.056 

-17.808 

-18.187 

-18.906 

-19.292 

-19.970 

-20.574 

-21.320 

17500 

57.234 

59.269 

61.502 

63.721 

65.983 

68.221 

70.482 

72.730 


-15.995 

-17.761 

-18.257 

-18.870 

-19.257 

,-19.938 

1 

-20.544 

-21.291 














while decreasing In others such that the total water vapor burden In the 

atmosphere under consideration remains unaltered. 

Computations of correction to the EBT were made for six multiples of 

the water vapor profile, namely, 0.0, 0.25, 0.5, 1.5, 2.0 and 3.0 in 

addition to the standard profile Itself, for all values of and the 

altitude. Table 4.5 lists the values of EBT and corresponding corrections 

for all values of T and the altitude of 17,500 ft. Results for this 

8 

altitude and T > 300, 310 €md 320”K were chosen for graphical illustration 
s 

and are shown in Figure 4.6. Polynomial regression in these results between 

Ax (water vapor burden) and AT (correction) reveals that cubic expressions 

are required to adequately represent these relationships. Regression 

coefficients a , a., a. and a. were determined for all cases and a 
o 1 z j o 

are found to be insignificant again. Coefficients a^, a 2 and a^ are 

presented in Table 4.6 and are entered together for each value of T^ 

and the altitude. It can be seen from this table that all these coefficients 

are functions of T as well as the altitude. Figure 4.7 shows the 
s 

variation of a^^, a^ and a^ with 7 for the altitude of 17,500 ft. 

Variation of these coefficients wi i altitude for T » 300*K is shown in 

s 

Figure 4.8. 

The effect of redistribution of water vapor in the atmosphere (without 
changing the total water vapor burden) has been studied by generating two 
modifications of the standard water vapor profile. Modification A is 
generated by arbitrarily decreasing the water vapor burden in the lower 
five layers and Increasing the burden by the same amount in the upper five 
layers. Modification B is generated by an exactly opposite operation. 

Table 4.7 shows a layer-by-layer distribution of water vapor burden for the 
standard and modified profiles. 
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Table 4.5 


Effective brightness temperature as a function of water 
vapor burden for an altitude of 17,500 ft. 


Actual 

Surface 

Effective Brightness Temperature (*K) 

Temp.(*K) 

0.00 

0.25 

0.50 

1.00 

1.50 

2.00 

3.00 

290 

290.000 

289.682 

289.348 

288.501 

287.404 

286.088 

283.060 


0.000 

-0.318 

-0.652 

-1.499 

-2.596 

-3.912 

-6.940 

295 

295.000 

294.561 

294.097 

292.910 

291.363 

289.523 

285.387 


0.000 

-0.439 

-0.903 

-2.090 

-3.637 

-5.477 

-9.613 

300 

300.000 

299.443 

298.856 

297.341 

295.354 

292.998 

287.759 


0.000 

-0.557 

-1.144 

-2.659 

-4.646 

-7.002 

-12.241 

305 

305.000 

304.330 

303.622 

301.786 

299.371 

296.515 

290.181 


0.000 

-0.670 

-1.378 

-3.214 

-5.629 

-8.485 

-14.819 

310 

310.000 

309.220 

308.394 

306.246 

303.414 

300.061 

292.645 


0.000 

-0.780 

-1.606 

-3.754 

-6.586 

-9.939 

-17.355 

315 

315.000 

314.113 

313.173 

310.719 

307.482 

303.638 

295.146 


0.000 

-0.887 

-1.827 

-4.281 

-7.518 

-11.362 

-19.856 

320 

320.000 

319.008 

317.957 

315.205 

311.566 

307.245 

297.688 


0.000 

-0.992 

-2.043 

-4.795 

-8.434 

-12.755 

-22.312 

325 

325.000 

323.908 

322.748 

319.700 

315.669 

310.875 

300.261 


0.000 

-1.092 

-2.252 

j 

-5.300 

-9.331 

-14.125 

-24.739 




Fig. 4.6 Variation of correction with water vapor burden from 17,500 ft. 
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Table 4.6 

Regression coefficients aj^, a 2 and a 3 as functions of actual 
surface temperature and altitude for water vapor burden study. 


Regression Coefficients 


290*K 295*K 300*K 305*K 310’K 315*K 320'K 325*K 


-0.025 

-0.015 


-0.088 

-0.054 

0.001 

-0.165 

-0.103 

0.004 

-0.327 

-0.215 

0.013 

-0.488 

-0.326 

0.023 

-0.635 

-0.410 

0.029 

-0.751 

-0.470 

0.034 

-0.841 

-0.511 

0.036 

-0.912 

-0.534 

0.037 

-0.983 

-0.551 

0.036 


-0.076 

-0.047 

0.000 

-0.216 

-0.133 

0.004 

-0.345 

-0.223 

0.010 

-0:579 

-0.392 

0.027 

-0.780 

-0.545 

0.044 

-0.943 

-0.664 

0.057 

-1.072 

-0.741 

0.065 

-1.173 

-0.790 

0.069 

-1.245 

-0.820 

0.071 

-1.322 

-0.839 

0.071 


-0.127 
-0.075 
.000 
-0.338 
- 0.211 
0.006 
-0.523 
-0.332 
0.014 
-0.820 
-0.563 
0.039 
-1.054 
-0.763 
0.064 
-1.241 
-0.908 
0.082 
-1.378 
-1.005 
0.094 
-1.486 
-1.063 
0.100 
-1.563 
-1.099 
0.103 
-1.641 
- 1.121 
0.104 


-0.173 
-0.105 
.000 
-0.457 
-0.283 
0.007 
-0.692 
-0.441 
0.018 
-1.047 
-0.735 
0.051 
-1.321 
-0.972 
0.082 
-1.524 
-1.149 
0.106 
-1.682 
-1.253 
0.119 
-1.797 
-1.321 
0.127 
-1.878 
-1.361 
0.131 
-1.960 
-1.386 
0.132 


-0.219 
.133 
.000 
-0.574 
.350 
.007 
-0.853 
-0.548 
0.022 
-1.270 
-0.898 
0.062 
-1.573 
-1.182 
0.100 
-1.807 
-1.375 
0.127 
-1.973 
-1.499 
0.143 
-2.094 
-1.576 
0.153 
-2.181 
-1.619 
0.158 
-2.265 
-1.649 
0.160 


-0.264 

-0.160 

0.000 

- 0.686 

-0.418 

0.008 

-1.008 

-0.654 

0.026 

-1.482 

-1.060 

0.072 

-1.827 

-1.376 

0.115 

-2.077 

-1.600 

0.147 

-2.254 

-1.739 

0.166 

-2.381 

-1.826 

0.178 

-2.473 

-1.872 

0.183 

-2.561 

-1.904 

0.186 


320'K 

325*K 

-0.307 

-0.349 

-0.187 

-0.211 

0.000 


-0.792 

-0.892 

-0.486 

-0.554 

0.009 

0.011 

-1.162 

-1.306 

-0.754 

-0.857 

0.029 

0.033 

-1.697 

-1.900 

-1.210 

-1.363 

0.080 

0.089 

-2.071 

-2.307 

-1.569 

-1.761 

0.129 

0.144 

-2.341 

-2.594 

-1.819 

-2.036 

0.166 

0.184 

-2.528 

-2.799 

-1.972 

-2.199 

0.188 

0.208 


-2.665 

-2.066 

0.200 

-2.763 

-2.114 

0.206 

-2.852 

-2.152 

0.210 


-2.942 

-2.301 

0.222 

-3.040 

-2.359 

0.229 

-3.134 

-2.397 

0.233 

















Effective brightnese tenperaturee have been obtained for all valuea 
of and the altitude of 17,500 ft. correaponding to the modified water 
vapor prof ilea and are presented along with the standard profile results 
in Table A. 8 The altitude of 17,500 ft. has been selMted here so that 
the observed effect is as large as possible. Comparison of the EBT values 
presented in Table 4.8 shows that except for the highest value of T^ , 
the differences between the EBT values corresponding to the standard and 
modified profiles are always smaller than 0.1*K. 


Table 4.8 

Effective brightness temperature for standard and 
modified water vapor profiles at 17,500 ft. 


Actual 

Surface 

Temp.CK) 

Eff. Brightness Temp. (*K) 

Standard 

Mod. A 

Mod. B 

290 

288.501 

288.453 

288.541 

295 

292.910 

292.883 

292.926 

300 

297.341 

297.334 

297.335 

305 

301.786 

301.798 

301.759 

310 

306.246 

306.277 

306.197 

315 

310.719 

310.768 

310.650 

320 

315.205 

315.271 

315.115 

325 

319.700 

319.784 

319.592 
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4.3 AtnoaphTic T«np«raturi Profllt 

Variations of atmospharic tanparaturt distribution ara axpactad to 
hava significant affact on tha upwalllng radlanca and tha EBT in this 
spactral ragion, particularly bacausa of tha strong tamparatura dapandanca 
of tha continuum absorption coafflciant [8,9]. This affact was axamlnad 
quantitatively by adopting four additional tamparatura profiles for the 
atmosphere. Since tha tamparatura gradient in tha troposphere is ralatlvaly 
I constant, profiles with fixed biases of 4-2, 41, -1, and -2*K ralatlva to 

tha standard wars considered as realistic variations. 

Effective brightness tamparatura for all values of and altitudes 
wars obtained corresponding to tha above taaiparature profiles. Since vari- 
ations of temperature profile become meaningless in absence of any water 
vapor, these computations have been carried out in the presence of a standard 
water vapor burden. EBT values for an altituue of 17,500 ft. are presented 
in Table 4.9. Corrections in this case are obtained relative to the EBT 
values which correspond to the standard temperature profile and are listed 
Just below the EBT values in the same table. The differences between T^ 
and T^ values which correspond to the standard temperature profile arc 
again due to the presence of water vapor. 

Figure 4.9 shows the correction as a function of the temperature bias for 

I two altitudes (4,500 and 17,500 ft.) and T ■ 300*K. Regression in these 

s 

data shows that these results can be adequately represented by linear 
expressions. Coefficients of regression a^ and a^ were determined for 
all cases. It is again clear from Fig. 4.9 that all these lines pass through 
^ the (o,o) point aiwl, therefore, the coefficients a^ arc insignificant. The j 

coefficient a^ for all values of and altitude are presented in Table 4.10. 

Figure 4.10 shows the coefficient a^ as a function of for the altitudes 

of 4,500 and 17,500 ft. It can be seen from this figure as well as from 



32 


i 


Tabl« 4.9 

Efftctivt brightntaa traparatur* for difftrtnt 
acaoapharic tanparatura proflla for cha altituda 
17,500 ft. 


Actual 

Surfaca 

E 

ffactiva Brlghtnaaa Taaparatura (*K) 

Taap(*K) 

AtBoapharlc Tamparatura Profilaa 



T' - T-2 

T' - T-1 

T' ■ T 


T* - T+2 

290 

288.235 

288.369 

288.501 

288.630 

288.758 


-0.266 

-0.132 


0.129 

0.257 

295 

292.644 

292.777 

292.910 

293.041 

293.169 


-0.266 

-0.133 

0.000 

0.131 

0.259 

300 

297.068 

297.206 

297.341 

297.474 

297.603 


-0.273 

-0..135 


0.133 

0.262 

305 

301.509 

301.649 

301.786 

301.920 

302.051 


-0.277 

-0.137 

0.000 

0.134 

0.265 

310 

305.965 

306.107 

306.246 

306.382 

306.515 


-0.281 

-0.139 


0.136 

0.269 

315 

310.433 

310.578 

310.719 

310.857 

310.992 


-0.286 

-0.141 

0.000 

0.138 

0.273 

320 

314.914 

315.061 

315.205 

315.344 

315.481 


-0.291 

-0.144 

0.000 

0.139 

0.276 

325 

319.405 

319.555 

319.700 

319.842 

319.981 


-0.295 

-0.145 

0.000 

0.142 

0.281 
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Table 4.10 


Regression coefficient as a function of actual surface 
cemperature and altitude for the temperature profile 
variation study. 


Altitude 

(ft) 

Regression Coefficients 

290"K 

295*K 

300*K 

305 *K 

310*K 

315 *K 

320"K 

325*K 

Average 

500 

0.017 

0.017 

0.016 

0.016 

0.015 

0.015 

0.015 

0.015 

0.016 

1500 

0.042 

0.042 

0.041 

0.041 

0.041 

0.041 

0.040 

0.040 

G.041 

2500 

0.062 

0.061 

0.061 

0.061 

0.061 

0.061 

0.061 

0.061 

0.061 

4500 

0.088 

0.088 

0.089 

0.089 

0.090 

0.090 

0.091 

0.092 

0.090 

6500 

0.106 

0.107 

0.107 

0.107 

0.109 

0.110 

0.111 

0.113 

0.109 

8500 

0.117 

0.118 

0.118 

0.120 

0.122 

0.123 

0.125 

0.126 

0.121 

10500 

0.123 

0.125 

0.125 

0.127 

0.129 

0.130 

0.132 

0.134 

0.128 

12500 

0.127 

0.128 

0.129 

0.132 

0.134 

0.135 

0.137 

0.139 

0.133 

14500 

0.129 

0.131 

0.132 

0.134 

0.136 

0.138 

0.140 

0.142 

0.135 

17500 

0.131 

0.131 

0.134 

0.135 

0.137 

0.140 

0.142 

0.144 

0.137 


Table 4.10 that the temperature dependence of the coefficient a^^ In this case 

is very weak. Average values of aj^ (over all values of T^) are obtained for 

each altitude. These are listed in the last column of Table 4.10 and shown In 

Figure 4.11. The figure shows a^^ as a strong function of altitude. 

The total correction in this case also consists of an additional component 

due to water vapor and can be expressed as 

AT * AT^ + AT (4.2) 

t w 

where AT^ is computed from the coefficients given in Table 4.10 and AT^ 
corresponds to the presence of standard water vapor burden. 
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5. MODEL CALCULATIONS FOR TWO PARAMETERS 


Deviation of only one parameter from Its initial value, as considered 
in the previous section, represents a highly idealized situation for the real 
atnwsphere where more than one parameter may undergo deviations at the same time. 
In this study we have examined some cases where two of the input parameters 
undergo deviations simultaneously. It is evident from the discussion in 
Section 4 that the Important parameters affecting the effective brightness 
temperature are: surface emlttance, water vapor burden, and the atmospheric 
temperature profile. Several combinations of these parameters have been examined 
in this section. It is Intended to establish relationships between the combined 
correction AT and the individual corrections AT caused by separate 

C X • 

variations of the two parameters. 

5.1 Surface Rmittance - Temperature Profile Variations 

Extensive model calculations were made to examine the relationships between 

the individual and combined corrections in this case. Computations were made 

for five values of surface emlttance, namely, 1.00, 0.95, 0.90, 0.85, and 0.80 

and five atmospheric temperature profiles. These temperature profiles were 

characterized by fixed biases of +2®K, +1®K, -1®K and -2**K in addition to the 

standard profile as given in Table 3.1. Computations were made for all eight 

values of T and calculations for an altitude of 17,500 ft. and T “ 325 ®K 
3 s 

are presented in Table 5.1. This particular set of parameters has been chosen 

for illustration to observe the largest differences between the sum of the 

individual corrections (Z AT ) and the combined correction AT . The first 

X c 

column of the table lists the values of surface emlttance and the deviation 
of emlttance from the initial value of unity. The next five columns list the 
values of the EBT and the differences (E AT^ - AT^) for the various temperature 
profiles. It can be seen from this table that the largest difference is 0.016®K 
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which Is much smaller than the acceptable 0.1 ‘K. It can be Inferred from 
these results that the combined correction In this case may be represented 
as the algebraic sum of the Individual corrections. Again, because of the 
presence of the standard burden of water vapor In the atmosphere In this case 
(without which the variations of temperature profile will be meaningless). It 
becomes essentially a three-parameter case. The total correction may be 
expressed mathematically as 

AT - AT + AT + AT (5.1) 

c X y w 

where ^T^ represents the correction due to standard water vapor burden 

and At & At represent corrections due to the deviations of the other two 
X y 

parameters. Also, It Is Important to note that correction AT^ due deviation 
of emlttance has to be calculated using regression coefficients obtained for 
the wet atmosphere (Table 4.4). 

The behavior of the corrections outlined abovu may be understood qualita- 
tively In physical terms as follows: the decrease of surface emlttance lowers 
surface radiance while the change in the atmospheric temperature profile 
affects the atmospheric radiance completely Independently. The change In 
absorption of the surface radiation by atmospheric constituents under the 
modified temperature profile has very small effect on the total upwelllng 
radiance and hence on the effective brightness temperature. 

5.2 Surface Emlttance - Water Vapor Burden Variations 

Model calculations In this case were performed for the same five values 

of surface emlttance as in subsection 5.1 and water vapor burdens of 0.0, 

0.5, 1.0 and 2.0. Results of this calculation for the altitude of 17,500 ft. 

and Tg * 325*K were chosen again for illustration and are presented in 

Table 5.2. The first column of this table lists the water vapor burden and 

the next five columns list the EBT and the differences (Z AT - AT ) for the 

X c 
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Table 5.2 


Slnultaneous Variation of Parameters 

Surface Emittance - Water Vapor Burden 

Altitude - 17,500 ft.; - 325 *K 

s 




various values of surface emittance. It is clear from Table 5.2 that the 



combined correction In this case cannot be represented by algebraic sum of 
the Individual corrections as In equation (5.1). Several analytical forms were 
tried for representing the combined corrections In this case and It was found 
after considerable numerical experimentation that It can be represented as 


AT - AT + AT + k, *AT *AT (5.2) 

c X y 1 X y 

the last term on the right-hand side representing the difference or residual 
correction. Values of the coefficient obtained for the various combinations 

of parameters Involved In Table 5.2 are listed just below the residual correc- 
tions. The blank spaces appear at those locations where the residual correc- 
tion as well as the Individual correction due to one of the parameters are 
both zero. For these cases, the coefficient k^ Is mathematically Indetermi- 
nate although practically zero. It can be seen from Table 5.2 that the value 
of kj^ for this set Is approximately constant. A mean value of 0.0197 + 0.003 
was obtained for this set. Similar mean values were obtained for all eight 
values of T^ and all ten altitudes, and are presented In Table 5.3. Figure 5.1 

shows k, as a function of T for z ■ 17,500 ft. (circles and solid lines) 

X s 

and also as a function of altitude for T ■ 300°K (squares and dashed lines) . 

s 

Value of k^ for any values of the parameters other than those listed In 
Table 5.3 may be obtained by linear Interpolation. 

It Is Important to note here that because water vapor burden Is Itself 
an Independent variable In this case and Individual correction due to water 
I vapor, AT Is calculated separately, the emittance correction, AT In this 

I ^ * 

case has to be computed using regression coefficients obtained for the dry 
^ atmosphere (Table 4.2). 



I 




Table 5.3 


Coefficient kj, representing residual correction for 
Surface Emlttance - Water Vapor Burden variations. 




ACTUAL SURFACE TEMPERATURE 



and with altitude (dashed curve) . 






The ohicrvtd differences or reslduel corrections in chl» cess oey be 
understood quelltetlvely in the folloving vey. The etaosphere being cooler 
then the surface for the cases considered, its effect is to attenuate the 
surface radiation. Since a decrease of jurface eaittance decreases the 
input surface radiance to the atmosphere, the net attenuation will be 
saaller in the latter case. This in effect is responsible for the over- 
prediction of correction observed in Table 5.2 when both deviations are 
considered separately. 

5.3 Water Vapor Burden - Temperature Profile Variations 

Model calculations in this case were carried out for the five tenpera- 
ture profiles mentioned in subsection 5.1 and three values of water vapor 
burden, namely 0.5, 1.0 and 2.0. Zero water vapor burden was not considered 
in this case again, because temperature profile variations loi>e meaning in 
absence of any water vapor. Results obtained for the altitude of 17,500 ft. 
and T^ " 325*R are presented lit Table 5.4 as an illustration, though they 
were obtained for all eight values of T^ and all ten altitudes. The first 
column of the table lists the water vapor burden (and the deviation of the 
water vapor burden relative to the standard as explained later) . The next 
five columns list the EBT values obtained for the various temperature profiles 
as indicated. It can be seen from Table 5.4 that for this case also the 
combined correction cannot be represented by the algebraic sum of the indivi- 
dual co-rectlons. The numbers entered below the EBT values are the residual 
corrections or the differences (I AT - AT ) . These differences are zero 

X c 

for T' • T because this represents essentially a one-parameter situation. 

The differences also equal zero for the standard water vapor burden because 
the one-parameter temperature profile variation calculations (subsection 4.3) 
were made in the presence of standard water vapor burden. Again, because 
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ch« ccoparatur* profile variations bacoaa naaninglass without watar vapor* 
standard watar vapor burdan is considarad as ona of tha naw initial conditions 
in this casa. Tha nuabars antarad balow tha watar vapor burden values in 
tha first column reprasant tha deviations of water vapor burdan relative to 
tha naw Initial value. 

Several analytical raprasantations ware tried and it was found that tha 
total correction for this sat of paraaatars may be rapresantad as 

AT - AT + AT + k, • AT' • AT (5.3) 

c X y 2 X y 

where AT and AT reprasant tha individual corrections due to watar 
X y 

vapor burdan and tamparatura profile deviations raspactlvaly. Tha last 
term again* raprasants tha residual correction and AT^ is defined as 

AX' - AT - AT (5.4) 

X w X 

where AT^ represents tha magnitude of correction due to standard watar 
vapor burdan. Tha values of the coefficient for tha cases covered in 

Table 5.4 arc entered Just below tha residual corrections in tha same table. 
Blank spaces occur at those places in this table where the residual correc- 
tion and ona of the individual corrections are both zero and* therefore, k^ 
is again mathematically indeterminate chough practically zero. Table 5.4 
also shows that k^ values for Ax ■ 2.0 (third row) are slightly higher 
than those for Ax * 0.5 (first row) while within each group the spread is 
very small. An average value of k^ " 0.249 0.014 was obtained for this 

set and it was determined that the maximum error introduced in the computation 
of AT^ using this average would still be much smaller than 0.1*K. Similar 
average values of k^ were obtained for all eight values of T^ and all ten 
altitudes and Che results are presented in Table 5.5. 

The magnitudes of the residual corrections as shown in Table 5.4 can also 
be explained qualitatively in terms of changes in attenuation of surface 
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Tab It 5.4 



SlBulttntout VtriaCion of Paraaatars 
Vatar Vapor Burdan - Tcmparatura Prof 11a 


Alt. - 17,500 ft.; T - 325*K 


Watar Vapor 
Burdan 


Ataoapharic Tamparatura Proflla 


T* ■ T-2 


T* - T-1 


T' 


T* ■ T+li T' - T+2 


0.5 

322.666 

i 

322.707 

322.748 

322.788 

322.827 

(-0.5) 

-0.213 

-0.104 

0.000 

0.102 

0.202 


0.237 

0.235 

- 

0.236 

0.236 

1.0 

319.405 

319.555 

319.700 

319. j;2 

319.981 

(0.0) 

0.000 

0.000 

0.000 

0.00;; 

0.000 

2.0 

309.893 

310.391 

310.875 

311.3 *5 

311.801 

(1.0) 

0.687 

0.339 

0.000 

-0.328 

-0.645 


0.264 


0.265 


0.262 


0.260 










ALTITUDE, Z 







Table 5.5 

Coefficient IC 2 representing residual correction for 
Water Vapor Burden - Temperature Profile variations. 


Altitude 

(ft) 

Actual Surface Temperature ("K) 

290 

295 

300 

305 

310 

315 

320 

325 

500 

25.260 

8.794 

5.540 

4.167 

3.489 

2.960 

2.664 

2.442 

1500 

7.319 

3.122 

2.079 

1.597 

1.320 

1.142 

1.028 

0.9353 

2500 

3.912 

1.947 

1.336 

1.049 

0.8743 

0.7642 

0.6832 

0.6226 

4500 

1.988 

1.159 

0.8461 

0.6782 

0.5753 

0.5056 

0.4522 

0.4124 

6500 

1.343 

0.8648 

0.6529 

0.5368 

0.4552 

0.4030 

0.3626 

0.3322 

8500 

1.063 

0.7208 

0.5587 

0.4606 

0.3987 

0.3533 

0.3190 

0.2924 

10500 

0.9199 

0.6448 

0.5098 

0.4232 

0.3680 

0.3278 

0.2963 

0.2722 

12500 

0.8383 

0.6002 

0.4785 

0.4018 

0.3499 

0.3126 

0.2827 

0.2588 

14500 

0.7886 

0.5774 

0.4596 

0.3983 

0.3397 

1 

0.3042 

0.2757 

0.2541 

17500 

0.7510 

0.5547 

0.4454 

0.3795 

0.3321 

0.2970 

0.2707 

0.2493 


radiation caused by changed water vapor burden and due to changes in 
emission from atmospheric water vapor because of the bias in the atmospheric 
temperature profile. 
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6. APPLICATION OF ALGORITHM 


The Information and algorithm developed In Sections 2, 4 and 5 can now 
be used to determine the correction to the measured effective brightness 
temperature for a known set of surface and atmospheric parameters. It can 
be evaluated using appropriate forms of equations (2.5) for one-parameter 
cases or (2.5) and (2.6) for more than one parameters. Various explicit 
forms of equation (2.5) were developed for specific cases as (4.1) and (4.2) 
and equation (2.6) as (5.1), (5.2) and (5.3). The regression coefficients 
required for these computations were determined and tabulated In Sections 4 
and 5. 


6.1 Evaluation of Correction 

It Is evident from Tables 4.2, 4.4, 4.6 and 4.10 that the regression 
coefficients a^ determined for various cases are dependent on the altitude 
as well as the actual surface temperature. Coefficients and k^ 

which are used In evaluation of the residual correction, listed In Tables 5.3 
and 5.5, also exhibit similar characteristics. Figures 4.5, 4.8, 4.11, 5.1 
and 5.2 show that all regression coefficients for any altitude between 500 
and 17,500 ft. can be obtained with desired accuracy by linear Interpolation 
from the tabulated values. Since the actual surface temperature, T^ Is the 
desired end result of the entire computation, an Iterative procedure has to 
be developed to arrive at its appropriate value. 

Values of the regression coefficients a as well as k. and k« for 

n i. z 

the appropriate altitude referring to T^ ■ 300“K are chosen as their initial 
estimates. Correction is computed using appropriate forms of equations (2.5) 
and (2.6). This correction combined with the measured value of EBT yields 
the first estimate of the actual surface temperature, say Values of 

the temperature-dependent coefficients a^ as well as k^^ and k2 are now 
obtained corresponding to found possible to represent che 
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temperature-dependence of each of these coefficients by analytical relations 
and, therefore, it was decided to obtain the coefficients corresponding to 
Tg^ by linear interpolation between tabulated values. Correction is evalu- 
ated again using these coefficients and equations (2.5) and (2.6). The 
second estimate of actual surface temperature, Tg^ is obtained by combining 
this latter correction with the measured EBT. This Iterative process is 
continued till successive estimates of Tg become consistent within 0.01”K. 

A computer program named CORRECT has been developed to accomplish the 

evaluation of correction using the above procedure. The program generates 

a table of regression coefficients for the altitude under consideration using 

interpolation, if necessary. A set of five parameters, namely, EBT (effective 

brightness temperature), ALT (altitude), EMI (surface emlttance), DW (water 

vapor burden) and DT (temperature profile bias) is supplied to this program. 

Depending upon the values of the last three parameters, the program decides 

if it is a one- or two-parameter case and automatically selects the appropriate 

equations for evaluating the correction. The program performs up to a maximum 

of ten Iterxatlons and if the T value does not converge as desired, it 

s 

reports the last estimate of T^ and indicates the nonconvergence. A listing 
of this program is also reproduced in Appendix B. 

6.2 Sample Calculations 

Extensive sample calculations were made to check out the working of 
program CORRECT using the EBT values obtained for various sets of parameters 
from the program PRTFIVE. Table 6.1 shows a set of such results covering a 
wide range of surface temperatures and other parameters. The first column of 
this table lists the values of the actual surface temperature (TEMS) used 
in the program PRTFIVE to obtaine EBT. The last column lists the values of 


the retrieved surface temperature (SURTEMP) obtained by the program CORRECT. 
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The other symbols used In this table have been explained earlier. The 
sets of parameters chosen for these calculations were intended to cover 
all possible combinations and also to show largest possible differences 
between the actual and retrieved surface temperatures. It is important 
to note that the combinations (1) surface emlttance - water vapor burden, 
and (il) water vapor burden-temperature profile (subsections 5.2 and 5.3) 
which Involve the computation of residual correction are likely to exhibit 
larger differences. A comparison of TEMS and SURTEMP values in this table 
shows that the largest difference is only 0.14‘R. 


TEMSCK ) 

EBT (® K ) 

ALT ( FT ) 

EMI 

DW 

DT (* K ) 

SURTEMP (* K ) 

325 

307.78 

10500 

0.80 

1 

0.00 

0.00 

325.00 

315 

296.48 

12500 

0.80 

1.00 

0.00 

314.99 

310 

297.33 

6500 

1.00 

3.00 

0.00 

309.97 

295 

293.21 

8500 

1.00 

1.00 

- 2.00 

295.00 

320 

315.60 

14500 

1.00 

1.00 

+ 2.00 

319.99 

305 

288.02 

17500 

0.80 

1.00 

- 2.00 

305.00 

310 

292.53 

17500 

0.80 

1.00 

+ 2.00 

309.97 

300 

282.88 

17500 

0.80 

2.00 

0.00 

300.10 

295 

294.00 

17500 

1.00 

0.50 

- 2.00 

294.94 

315 

304.51 

17500 

1.00 

2.00 

+ 2.00 

314.86 

320 

295.36 

17500 

0.80 

2.00 

0.00 

320.04 

315 

312.26 

17500 

1.00 

0.50 

+ 2.00 

314.92 


















7. SENSITIVITY CALCULATIONS 


It Is important in a study of this type to examine the sensitivity 
of the retrieved surface temperature to the uncertainties in the various 
input parameters. It is possible In cases of certain parameters to 
estimate the 'sensitivity from the information already presented in Section 4. 
For others, however, separate calculations had to be made for this purpose. 
Effects of uncertainty in the values of surface emlttance, water vapor 
burden, temperature distribution and altitude of observation are examined 
in this section. 

7.1 Surface Emlttance 

Sensitivity of the retrieved surface temperature to uncertainties in 
assumed value of surface emlttance may be estimated from the results pre- 
sented in Tables 4.1 and 4.3 for dry and wet atmosphere conditions, respec- 
tively. These tables list the values of EBT as a function of T^ for several 
values of surface emlttance. Differences between EBT values referring to 
e * 1.00 and 0.95 are indicative of the differences between the measured 
values of EBT expected as a result of a five percent variation of surface 
emlttance. Since surface temperature is obtained directly from the EBT, it 
is reasonable to assume that these differences are approximately equal to 
the errors Introduced in the retrieved surface temperature due to five percent 
uncertainty in the assumed value of surface emlttance. 

Figure 7.1 shows EBT as a function of the actual surface temperature 
obtained with two values of input emlttance (1.00 and 0.95). The dry atmo- 
sphere results are Independent of altitude and the wet atmosphere results 
refer to the altitude of 17,500 ft. It is evident from figure 7.1 that for 
both cases (dry and wet) the differences between the two EBT curves are 
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ATMOSPHERE 


1.00 


0.95 


/ 

/ 


/ y / 
/ / / 

# 


WET 

ATMOSPHERE 
ALT = 1 7500 FT 


>0 300 310 320 325 

ACTUAL SURFACE TEMPERATURE (®K) 

Fig. 7.1 Data showing sensitivity of EBT to uncertainty in the 
value of surface emittance. 





substantial (an average of approx. 3.3*K). It is Inferred, therefore, that 
surface emlttance is a very important parameter and great caution should 
be exercised in assuming its value in the present work. 

7.2 Water Vapor Burden 

Sensitivity of the retrieved surface temperature to uncertainties in 
the assumed water vapor burden has been estimated by computing EBT for 
water vapor burdens of 0.75 and 1.25 and comparing the results with EBT 
values obtained for standard water vapor burden. Table 7.1 shows results for 
all eight values of T and the altitude of 17,500 ft. The results presented 

3 

in this table are depicted in Figure 7.2 and show that a 25 percent varia- 
tion of water vapor burden affects the measured EBT from 0.5 to 1.5*K 

(average of about 1*K) for the range of T covered in the present work. 

8 

It is easy to visualize that the values of surface temperature obtained from 
these values of EBT will differ by an average of approximately 1*R. It can 
be Inferred, therefore, that an uncertainty of 25 percent in the waver vapor 
burden Introduces an error of approximately I’K in the retrieved surface 
temperature. 

* 

7.3 Atmospheric Temperature Distribution 

It can be seen from Table 4.9 that a constant bias of 2*K on the atmo- 
spheric temperature profile affects the EBT measured from 17,500 ft. by an 

average of approximately 0.27®K for the range of T under consideration. 

8 

As discussed in the previous subsections, it can be inferred from these 
results that an uncertainty of 2®K in the atmospheric temperature distribution 
throughout the atmosphere introduces an error of approximately 0.27*K in the 
retrieved surface temperature. 
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EFFECTIVE BRIGHTNESS TEMPERATURE 






Table 7.1 


Results showing the effect of uncertainty in 
water vapor burden on effective brightness 
temperature. 


Actual 
Surface 
Temp. (*K) 

Effective Brightness Temp. (*K) 

Water Vapor Burden 

JT73 


1.25 

290 

288.96 

288.50 

287.98 

295 

293.55 

292.91 

292.18 

300 

298.16 

297.34 

296.40 

305 

302.78 

301.79 

300.65 

310 

307.41 

306.25 

304.91 

315 

312.05 

310.72 

309.19 

320 

316.70 

315.21 

313.49 

325 

321.35 

319.70 

317.80 


7.4 Altitude of Observation 

rhe effect of uncertainty in the altitude of observation on retrieved 
surface temperature is estimated by computing EBT as observed from 17,000 
and 18,000 ft. and comparing the results with those obtained for 17,500 ft. 
The differences observed for both cases were very small (£0.02*K) and were 
not considered significant. It was Inferred, therefore, that an uncertainty 
of up to 500 ft. (at an average altitude of 17,500 ft.) In the altitude of 
observation introduces no significant error in the retrieved surface 
temperature. 
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8. CONCLUSIONS 


Th« praatnt study shows that It is posslbls to dstsrnias ths stao- 
sphsric corrsctlon to ths sffsetlvs brightnsss tsapsrsturs ss nsssursd 
rsDOtsly by s rsdistlon thernoiBStar without having to psrfora dstallsd 
radlatlvs transfsr calculations. It Is shown from ths nodal calculations 
in Section 4 that corrections can be represented accurately by simple 
analytical expressions (linear, quadratic or cubic) in terns of ths 
magnitudes or deviations of the corresponding surface or atmospheric 
p.'irameters. It is also shown from the model calculations in Section 5 
that combined correction for simultaneous variations of several paranwters 
can be represented in terms of their Individual corrections. It has been 
shown further (Section 6) that these simple expressions can be used to 
compute the corrections with a high degree of accuracy from a knowledge of 
the magnitudes or deviations of the various surface end atmospheric para- 
meters over wide range conditions. The largest cumulative mathematical 
error introduced by this algorithm was found to be less than 0.1S*K. It 
can also be seen from the results of the sensitivity calculations presented 
in Section 7 that uncertainties in the asswaed values of surface emittance 
and water vapor burden introduce large errors in the retrieved surface 
temperature. Uncertainties in the atmospheric temperature distribution and 
the altitude of observation are found to have relatively smaller effects. 

It is estimated that the use of this algorithm instead of the expenalve 
radiative transfer calculations will reduce the time and cost Involved in 
surface temperature retrlev.:! by at least an order of magnitude. Although 
the ^results presented here refer to one particular instrument (PRT-S) , the 
technique is general and should be applicable to any radiometric measurement. 
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APPENDIX A1 




A 

AC 

ACB, ACE 
& ACM 

AL 

ALFT, 

ALKM 

ALX 

ATR 

B 

BETA 

BTR 

DEL 

DLIM 

EL 

EMI 

FACT 

FIL 

FR 

FRB 


EXPLANATION OF SYMBOLS USED IN 
PROGRAM "PRTFIVE" AND ITS SUBROUTINES 


Constant used In representing frequency variation of the 
continuiim absorption coefficient. 

Total monochromatic absorption coefficient. 

Wing contribution to the total absorption coefficient at the 
interval boundaries and center. 

Array of half-widths of the individual lines. 

Altitude of observation in feet and kilometers, 
respectively. 

Width of a narrow sub-interval near line center. 

Integrated upwelling radiance at the altitude of observation 
including atmospheric attenuation. 

Another constant used in representing frequency variation 
of the continuum absorption coefficient. 

Exponent used for the above purpose. 

Integrated upwelling radiance at the altitude of observation 
excluding atmospheric attenuation. 

w!dtb of an interval. 

Width of the region from which the direct contribution 
is obtained. 

Array of energies of the lower states for individual lines. 
Surface emlttance. 

Factor used for obtaining line-intensities corresponding 
to the temperature of each layer. 

Array of filter function for all Intervals. 

Array of frequencies of the individual lines. 

Frequencies at the boundaries of intervals. 


FRC 

FRE 

FSG 

FRS 

FRT 

IG 

LA 

L£ 

LLR 

LX 

MP 

NQ 

NT 

PART 

PL 

PLK 

PNTP 

PPL 

PREC 

PSK 

PW 

PXK 

QV 

RAD 


Frequencies at the centers of intervals. 

Frequencies of all lines within one Interval. 

Frequencies of all locations where absorption coefficient is 
calculated within one Interval. 

Frequencies at the sub-interval boundaries within one Interval. 
Tabulated values of filter function. 

Total nunber of frequency locations at which transmittance 
calculation is made, within one interval. 

Number of layers for a given altitude of observation. 

Total number of lines in the entire frequency range. 

Number of different altitudes of observation. 

Array of n..jiber of layers for different altitudes of observation. 
Number of lines within one Interval. 

Number of sub-intervals in an Interval. 

Number of surface temperature values for which calculations 
are made. 

A lumped constant accounting for the variations of vibrational 
and rotational partition functions. 

Monochromatic optical thickness. 

Planck functions for different atmospheric layers. 

Pressure at NTP (760 mm Hg) . 

Pressure path length for each layer. 

Average pressure for each layer. 

Planck functions for dlffert t values of surface temperature. 
Partial pressure of water vapor for each layer. 

Planck functions for surface for the temperature array TEMX. 
Concent rat loo, of water vapor in different layers (in ppmV) . 
Upwelllng radiance in one interval at one particular temperature. 


63 


RP 


Exponent which accounts for the temperature dependence of 
the rotational partition function. 

Array of Intensities for Individual lines. 




SI 

TEMC 

TEMR 

TEMS 

TEMX 

TEMY 

THC 

TNTP 

TR 

TRA 

TRS 

VPF 

WF 

WIDF 


Average temperature for each layer. 

Reference temperature for the water vapor line parameters. 

Array of surface temperatures for which calculations are made. 
Array of surface temperatures which are used for Interpolation. 
Computed values of effective brightness temperatures. 

Thickness of each layer. 

Temperature at NTP (273*K). 

Monochromatic transmittance at each of the frequencies FRG. 

Averaged transmittance for each Interval for different layers. 

Averaged transmittance for each Interval In the Subroutine TRANS. 

Factor accounting for the temperature dependence of the vibra- 
tional partition function. 

Frequencies at which the filter function is tabulated. 

Factor used to convert line half-widths from reference values 
to those appropriate for each layer. 


WLl, WL2 

WLIM 

XL1.XL2 


Weighting factors used in the Gauss-Legendre quadrature formula. 
Width of the regions from which the wing contribution is obtained. 
Abscissa values used in the Gauss-Legendre quadrature formula. 






APPENDIX A2 


A 

A1.A2. 

& A3 

AAl, AA2, 
& AA3 

AD 

ALFT 

ALT 

ARI, AR2, 
& AR3 

ASTM 

ASTN 

AW, AX, 

& AY 

AZl, AZ2, 
& AZ3 

B 

C 

Cl, C2 

CZl, C22 

DE 

DT 


EXPLANATION OF SYMBOLS USED IN PROGRAM "CORRECT." 


Current values of the regression coefficients for the first 
parameter . 

Regression coefficients for all altitudes, temperatures and 
parameters. 

Regression coefficients for all altitudes and temperatures 
for one particular parameter. 

Regression coefficients for the dry atmosphere surface 
emittance variation case. 

Altitudes at layer boundaries In feet. 

Given altitude of observation. 

Regression coefficients for the given altitude of observation 
for all temperatures and one parameter. 

Initial estimate of the surface temperature (300°K In all cases) . 

Current estimate of surface temperature during Iteration. 

Regression coefficients for water vapor, first parameter and 
second parameter, respectively, for the given altitude and all 
temperatures. 

Regression coefficients for the given altitude, all temperatures 
and all parameters. 

Current values of the regression coefficients for the second 
parameter. 

Current values of the coefficients for the residual correction 
terms. 

Coefficients for the residual correction terms for all altitudes 
and temperatures. 

Coefficients for the residual correction terms for the altitude 
of observation. 

Deviation of surface emittance from unity. 

Bias In atmospheric temperature profile. 


Atmospheric water vapor burden in multiples of standard. 


DX, DY 


Deviations of the two parameters under consideration. 


Measured effective brightness temperature. 


Surface emlttance. 


Convergence parameter. 


SURTEMP 


Retrieved surface temperature. 


Array of surface temperatures for which calculations are made. 
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APPENDIX B1 


LISTING OF COMPUTER PROGRAM "PRTFIVE" 
AND ITS SUBROUTINES 


C 

C 


PROGRAM PRTF 1 VE U NPUT « OUTPUT . T APE2 ) 

INTEGER X«Y»T 

0 1 MENS I ON AlFT ( 1 0 > ♦ AlKM ( 1 0 ) « THC (10). PREC (10).LX(10)» 
/VPF(IO) tOVnO) .PPL! 10) «FRT(25) .K/F(25 ) »FGR<24 ) .FIL ( 120) . 
/ATR<8» 10 )»BTR( 1 0) .FRE( 15) .PL <200 ) . PSK < S ) . PL< < 1 0 ) . PXK ( 1 0 ) « 
✓TEMS(8) .TEMX (10 ).TRA (11) .RAO <8 ) .TEMY (8. lO ) . lER (8. 1 0 ) 
C0MM0N/A8S0RB/FR(320 ) .SI ( 320 ) .tL <320 ) .AL t 320 ) . I- R(i t 2U{) ) . 
/AC ( 200 ) . FR3 (121 ) . FRC ( I 20 ). W I OF ( 1 0 ) . FACT < I 0 ) . PAR t (10). 
/TEMC( 10 ) .PW < 1C ) .OLl V ,*L1M,0ELA.8 ETA.TE-:mR. TO.LF . IG.PI .A.B 
C0MM0N/TRANE/FRS(51 ) .TR(200 ) .WLl .WL2.NQ.DEL 


READ 

10. 

EMI 

READ 

10. 

FRl.FRU.DEL 

READ 

1 1 . 

LE.NT. JT.LLH 

READ 

10. 

dlim.wlim,alx 

READ 

10. 

PNTP . TNTP ♦ TEMR , TO 

read 

10. 

RP 

READ 

12. 

A,8.8ETA 

READ 

13. 

(AlFT(LL ) ,LL=1 .LLH) 

READ 

14. 

(PREC(L) «L=1. 10) 

READ 

10. 

(TEMC(L) ,L=i ,10) 

READ 

1 0 . 

( TEMSt T ) ,T=1 ,nT) 

READ 

10. 

(TEMX( J) ,U=1 , JT) 

READ 

14. 

(THC<L ) ,L=1 , 10 ) 

read 

15. 

(FRT (M ) ,M=1 ,25 ) 

read 

14, 

( VPF (L ) ,L= 1 ,10) 

read 

16. 

XLl ,XL2, WLI , v,l2 

READ 

12. 

(CV iL ) .L= 1 » I 0 j 

read 

(2.17) (FR(X ) .SI (X ) ,EL (X ) .AL (X ) .X = l ,LE ) 


10 FORMAT < 10F8.2 ) 

11 FORMAT! 1615) 

12 FORMAT (5E16. 3 ) 

13 FORMAT ( 1 0F8.0 ) 

1 4 FORMAT ( 1 0F8.4 ) 

15 FORMAT! 10F8.3) 

16 FORMAT (8^10.6 ) 

17 FORMAT(2(Fl0.3.cl2.4,F10.3.F5.3.3X) ) 

determines the number Cf- LAYtRb UW EACH U^- I me 

several altitudes Or OBSERVATIONS 
DO 1 00 LL=l .LLH 

AlKM <LL ) =ALFT (LL ) *3 . 04BE-04 +0 • U5 
L = 0 

STHC-0. 


101 L=L+1 

STHC=STHC + ThC (L ) 




;.vU 




Ui'.iUi N A.ij i’Abr. 

OF POOR QUALITY 


ALR« ALKM ( LL ) -STHC 
IF (ALR.LT.O) GO TO 100 
LX<LL)»L 
GO TO 101 
100 CONTINUE 

COMPUTES PARTIAL PRESSURE t PRESSURE PATH LENGTH 
FOR water vapor AND SOME OTHER ALTITUDE 
DEPENDENT PARAMETERS 
PI=3*14159 
AWRsl.6 

C0NST=1 .E+05*TNTP/PNTP 
DO 102 L=1 « 10 

PART ( L ) *VPF < L ) ♦ ( TEMR/TEMC ( L ) ) **RP 
FACTIL )=1 .439* (TEMC (L )-TEMR)/(TEMC{L )*TEMR) 
W10F(L)=(SQRT(TEMR/TEMC(L) ) )*PREC (L )/PNTP 
PPL(L)=CONST*QV<L )*PREC(L )*THC (L)/TEMC(L ) 

102 PW(L)=QV(L)*PREC(L) 

DIVIDES the entire FREQUENCY RANGE INTO NARROW INTERVALS 
DELA=0.5*DEL 
R<={FRU-FRL)/DEL+0, 1 
<RsRK 

FRB( 1 )=FRL 
DO 103 K=1 .KR 
FRC(K )=FRa(K )+DELA 

103 FRB(K+1 )=FRB (tC )+OEL 

COMPUTES THE INSTRUMENT FILTER FUNCTION FOR 
EVERY interval FROM THE TABULATED VALUES 
WF { 1 )=FRl 
DO 104 M=1.24 
WF(M+1 )=KF(M) + lO, 

104 FGR(M)= (FRT(M+1 )-rRT{M) )/(WF(M+l )-wF<M) ) 

DO 105 <=1 «<R 

M=0 

106 M=M+1 

IF (FRC (< ).LT.WF(M) ) GO TO 105 
IF (FRC(K ).GE.WF(M+1 ) ) GO TO 1 U6 
FlL(0 = rRT(y)+FGRfM)*( FRC { <) - WF { M ) ) 

105 CONTINUE 
ALY = 2.*Ai_X 

CONS= 1 6 » *6 .6255-0 7 
CNST=6. 625*0. 3/1 .38 

initializes total UPWELLING radiances for DIFFERENT 
ALTITUDES AND SURFACE TEMPERATURES 
DO 107 LL=l .LLH 


<^0 J07 T=1 »NT 

107 ATR(T»LL)»0. 

00 lOa 

108 BTR(J)»0. 

Y»0 

C STARTS COMPUTING TRANSMITTANCE AND UPWEULING 
C RADIANCE FOR ONE INTERVAL AT A TIME 
DO 109 K*1 »KR 
M=0 

C DETERMINES THE NUMBER OF LINES IN ONE INTERVAL 
MP=M 

111 IF lY.GE.LE) GO TO llO 
Y*Y+1 

IF <FR(Y) .LT.FRB(K) > GO TO 111 
IF <FR(Y).GE.FRB(K+1 M GO TO 112 
M=M+1 

FRE<M)=FR(Y) 

MPsM 

GO TO 1 11 

112 YsY-1 
110 CONTINUE 

N*1 

C determines THE NUMBER OF SUB- I NTERVAlS IN ONE INTERVAL 
FRS{N)*FR&{< ) 

IF (MP.LE.O) GO TO 113 
DO 1 lA M-*l .MP 
OIFaFRE «M )-FRS (Nl 
IF (DIF.LE.ALX) GO TO 115 
IF (DIF.LE.ALY) GO TO 116 
FRS(N+1 )=FR£(M)-ALY 
FRS(N+2 )=FRE(M)-ALX 
FRS(N+3)=FRE(M) 

N=N+3 
GO TO 117 

116 FRS(N+1 )=FRE(M)-ALX 
FRS(N+2 )=FRE(M ) 

NsN+2 

GO TO 117 

1 15 FRS(N+1 )=FRE(M) 

N = N+ I 

117 IF (M.EC..MP) GO TO 1 18 
OlF=FRE(M+l )-FRS(N) 

IF (OIF.LE.ALX) GO TO llA 
IF (01F.GT.ALY) GO TO 119 






FRS(N+l )=FRE(M)+ALX 
N«N4>1 
GO TO 114 

119 FRS(N+1 >»FRE(M)4-ALX 
FRS ( N+2 ) »FRE ( M )+ALY 
N«N-f2 

114 CONTINUE 

lia DlF*FRd(K+l 1 -FREcM) 

IF (DIF*LE»ALX) CO TO 113 
IF (DIF.LE.ALY) GO TO 120 
FRSIN+l )=FRE(M)+ALX 
FRS(N+2 )*FRE(M l+ALY 
FRS(N+3 )aFR8<K+l ) 

NPsN+3 
GO TO 121 

120 FRS(N+1 )»FRECM l+ALX 
FRSIN+2 )=FRB (K+1 ) 

NP*N+2 

GO TO 121 

113 FRSIN+l )=FRB (K+1 ) 

NP»N+1 

121 NQ=NP-1 

DETERMINES THE FREQUENCIES AT ALL GRID POINTS WlTHlN ONE INTERVAL 
IG=4*NQ 
DO 122 Nsi ,N0 
I a4*(N-l ) + l 

VARaO.5* (FRS(N+1 )-FRS(N) > 

CON=0,5*(FRS(N+1 )+FRS(N) ) 

FRGII )3C0N-V^R*X'_ 1 
FRG(1+1 ) =C0N-VAR*XL2 
FRG( 1+2 )»C0N+VAR*XL2 

122 FRG( I+3)»C0N+VAR*XL1 
PNUMaDEL*CONS*FRC (< )+FRC(K )*FRC <K )*1 .E-07 
EEX=CNST*FRC<K ) 

COMPUTES values of PLANCK FUNCTIONS FOR DIFFERENT VALUES OF 

SURFACE temperature AnO DIFFERENT LAYERS OF THE ATMOSPHERE 
DO 123 T=1 .NT 

123 PSK(T)sPNUM/(EXP(EEX/TEMS{T ) )-l . ) 

DO 124 L=1 . 10 

124 PLK(L )=PNUM/(Exr vl.EX/TEMC(L ))-!•) 

DO 125 J=1 . JT 

1 25 PX< ( J 5 =PNUM/ (EXP( EEX/TEMX ( J ) )-l . ) 

DO 132 J=l.JT 

132 8TR(J)=BTR( J)+PXK (J)*FIL(K) 




OF 
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n o 


C STARTS TO WORK FOR EACH ALTITUDE SEPARATELY 
DO 109 LL»1 «LLH 
DO 126 I>1*IG 

126 PL (I )=0. 

LA«LX<LL ) 

LB«LA-t-l 
TRA(LB)=1 • 

C CONSIDERS EACH LAYER SEPARATELY. CALLS SUBROUTINE 
C COEFF TO evaluate ABSORPTION COEFFICIENT AT EACH 
C GRID POINT. Then computes the optical DEPTH AND 

c monochromatic transmittance 
DO 127 M=1.lA 
LsLA+I -M 
call COEFF(L.K) 

DO 128 I >1 . IG 

PL ( I ) =PL < I ) +AC ( I ) *PPL ( L ) 

IF (PL ( I ) .GT.675. ) GO TO 129 
TR( I )»EXP(-PL( I ) ) 

GO TO 128 

129 TR( I )»0. 

128 CONTINUE 

CALLS subroutine TRANS TO COMPUTE AVERAGE 
TRANSMITTANCE FOR EACH INTERVAL 
call TRANS(TRA(L) ) 

127 CONTINUE 

c COMPUTES Integrated upwelling radiance for 
c different altitudes and Surface temperatures 
COMP=0* 

DO 130 L»1.LA 

130 COMP*COMP+PLK '•_)»( TSA (L+1 '-TRAIL: ) 

DO 131 T*1.NT 

RAD(TJ=C0MP+EMI*DS<(T)*TRA( 1 > 

131 ATR(T.LL)=ATR(T.LL>+RAD(T)-fFIL(K) 

109 CONTINUE 

NMAX=1 0 
NN=1 0 
NTAB= 1 
I0R = 2 

C USES THE INTERPOLATION ROUTINE lUNl TO EVALUATE 
C THE EFFECTIVE BRIGHTNESS TEMPERATURE 
DO 133 LL=1 .LLH 
DO 133 T=1 .NT 
XO=ATR ( T .LL ) 

IPTr- 1 




CALL lONl (NMAX*NNfBTR»NTAB*TEMX» lOR*XO«YO« IPT * lERR) 

TEMY(T«LL)»VO 

lER(TtLL)»lERR 

133 CONTINUE 

DO 134 LL«1«LLH 

PRINT 60. ( ',‘EMS(T ) .ATRIT.LL ) »TEMY<T.LL J . lER(T.LL) »T«1 .NT ) 

134 CONTINUE 

PRINT 61 » (TEMX(U).BTR(J)»J»1»JT> 

60 FORMAT! IHl//// (Fl5*2 .E20.5»F15*3. 110.//) ) 
e»l FORMAT! 1H1////<F15.2«E20. 5*//) ) 

STOP 

END 






0 ? 


72 


o o 


subroutine COEFF(U.K) 

INTEGER X 

DIMENSION CONT(EOO) 

C0MM0N/ABS0RB/FR(320 )*SI ( 32U ) « EL ( 320 ) • AL ( 320 ) « FRG ( 200 ) « 
/AC (200) .FRa( 121 ) »FRC< 120) . jk1DF( 10) .FACT ( 10) .PART ( 10 
/TEMC( 10 ) «PW( 10 > .DLIM.WLIM.OElA.OLTA.TEMR. TO .LE . lO.P I .A.B 
C INITIALIZES ALL ABSORPTION COEFFICIENTS TO ZtHO 
ACB>ACM>ACE>0. 

DO 300 1*1 . IG 

300 AC( I )*0. 

C STARTS TO GO THROUGH ThE LINES TO DETERMINE WHICH 
C ONES CONTRIBUTE FOR THE PARTICULAR INTERVAL 
DO 301 X*1 .lE 
01F«ABS(FR(X)-FRC(<) ) 

IF (OIF.GT.WLIM) GO TO 301 

S1A«SI (X)*PART(L)*EXP(EL(X )*FACT (L) ) 

ALB»AL(X)*W1DF(L) 

IF (OIF.GT.DLIM) GO TO 302 
C COMPUTES THE DIRECT CONTRIBUTION 
DO 303 1*1 . IG 
FO«FR(X)-FRG( I ) 

DEN«PI * (FD*FD+ALB»AlB ) 

3C3 AC( I )*AC( 1 )+SI A*ALB/DEN 
GO TO 301 

evaluates the WING CONTRIBUTION AnQ ADOS 
THE DIRECT and WlNG CONT^V J 3UT I ONi 
302 F8*FR(X )-FRB (< ) 

FM»FR(X )-FRCt< ) 

FEaFR( X )-FRB (<+l ) 

PNUM*S I A*ALB 
ACB*ACB+PNUM/ (P I *F5»F5 ) 

ACM»ACM+PnUM/(PI*FM*FM ) 

ACE«ACE+PNOM/(PI*FE*='£ ) 

301 continue 

SLla{ACM-ACB)/(FRC(<)-FRB(K) ) 

SL2a<AC£-ACM)/(FR8<K+l )-FRC(< ) ) 

DO 304 1*1 . IG 
DIF»(FRG( I )-FRB<< ) ) 

IF (DIF.GE.DELA ) GO TO 3C5 
AC( I )»AC( I )-*-ACS + SLl *OIF 
GO TO 304 

305 AC( I )*AC< I ) +ACMfSL2* (DIF-DELA ) 

304 continue 

C COMPUTES THE CONTINUUM PART OF THE ABSORPTION COEFFICIENT 
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C AND ADOS TO THE LINE PART TO EVALUATE THE TOTAL 

FAC»EXP<T0*(TEM«-TEMC(L) )/<TEMC (L>*TEM«» >#P\k{L» 
00 306 1 ■ 1 * I G 

C0NT( 1 )b(A4>B»ExP(-B£TA*FRG( I ) ) )*FAC*2.69E + 19 
306 AC(I )-AC( n-»>CONT( I ) 

RETURN 

END 


subroutine TRANS (TRS) 

C0MM0N/TRANE/FRS(51 ) •TR(200) «wL1 «WL2*NQ«0EL 

c evaluates the average transmittance for each 
C SUB- interval using FOUR-POINT GAUSS-LCGENORE 
C QUADRATURE FORMULA ANO THEN COMPUTES THE 
C AVERAGE TRANSMITTANCE FOR THE ENTIRE INTERVAL 
TRS«0* 

DO 400 N>1 *N0 
I «4*<N-1 J4-1 

VAR«0.5*(FRS(N+1 )-FRS(N) > 

SUMt«TR( I )+TR( 1+3) 

SUM2«TR( I + I >+TR( 1+2 ) 

SUM-WLl *S'JM1 +4L2*SUM2 
400 TRS»TRS+SUM*vAR/OEL 
RETURN 
END 
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APPDfDIX B2 


LISTING OF COMPUTER PROGRAM "CORRECT” 


C 

c 


c 

c 


PROGRAM CORRECT < I NPUT » OUTPUT . TAPE2 ) 

DIMENSION ALFT( 10) *T£M5(a) * A1 U0«e«3)*A2U0«a*3)« 
/A3(iO«8*3)«Cl ( I0«d) tCEl 10*8) «AZ1 (8*3) «AZ2(8«3) 

/CZl (8) oCZP(8) tAAl (10,8) «AA2( I 0 1 8 ) , A A3 d 0 , 8 ) • AR 1 ( 8 ) • 

/AR2 (8) t AR3(8) « AX ( 8* 3 ) • AY ( 8 * 3 ) * AW ( 8 * 3 ) * AQ ( 8* 3 ) • A ( 3 ) « 

/B(3) «C(3) * 1ERL(3,3) * l£RC( IC ) * 1£RX( lO) « l£RY (10) , IERW( 10) 

DATA ICON* lERCl • 1ERC2* IERL*ASTM/1 ,l 1»9*300./ 

DATA lERC. !ERX« lERY* IERW/40*9/ 
data ALFT/SOC* * 1500* *2500. *4500* * 

/6500* «6500* , lOSOo* , 12500* , 14500* ,17500./ 

DATA TEMS/29C..29S. ,300* . 

/305* *310. ,315* ,320* ,325*/ 

READ (2* 10) (((A1(L*J*K)*L> 1*10)* JB 1,0), 

/( (A2(L*J*<)*L«1*10)*UB1,8>*((A3(L*J,K)*L*1*10)«J«1 *8),<>1 *3) 
READ (2* 10) ((Cl (L«J),L»1,10)*J»1.8),((C2(L*J>*U»1»10),J«1,8) 
READ (2* 11) ( (A0( J*N) , J-1 *8 ) ,M«1 • 3 > 
read 12, EoT.AlT.EMI ,D*,DT 

10 FORMAT ( 10P8.4 ) 

11 FORMAT (8F8.4 ) 

12 FORMAT ( 10F8.2 ) 

DE-EMl-1 . 

1 ORDER* 1 

determines if the altitude OF Observation is 

ON ONE OF The LAYER BOUNDARIES 
00 100 L-1 * 10 
ADF»0*C2*AlFT (L ) 

DlFF«AbS ( ALFT (L )-AlT ) 

if(DIff.gt.aof) go to luo 
GO TO 101 

100 CONTINUE 
GO TO 102 

IF ON THE layer boundary, USES THE COEFFICIENTS 

FOR that altitude as They are 

101 ILBL 

DO 103 K»l*3 
DO 103 J*1 *8 
AZl ( J,K )»Al ( IL, J*< ) 

AZ2( J*K)>A2( IL,J*<) 

AZ3( J*< )=A3( IL,Ut<) 

103 continue 

DO 104 J=l,8 
CZl (J)*C1(IL*0) 

CZ2( J)*C2( IL»U) 

104 continue 






c 

c 


c 

c 

c 


GO TO 139 
102 NMAX^NNalO 
NTAB=a 


ORIGINAL PAGE IS 
OF POOR QUALITY 


IF NOT* GENERATES TABLE OF COEFFICIENTS AT THE 
altitude OF OBSERVATION BY LINEAR INTERPOLATION 
00 105 K*1 .3 
DO 106 U*1 »8 
DO 106 L=1 tlO 
AAl (L* J)=A1 (L«J«K ) 

AA2(L» J)=A2(L«J»K ) 

AA3(L» J)=A3(L»J«K ) 

106 continue 


1PT1=IPT2=IPT3=-1 


CALL lUNI (NMAXtNN ALFT«NTAa» 
/AA 1 « 1 ORDER » ALT , AK 1 . I PT 1 ♦ lERR ) 
IERLU tK) = IEPR 

CALL lUNI (NMAX.NN.ALFT^NTAB. 
/AA2 , I ORDER . alt . AR2 « I PT2 » I ERR ) 
IERL(2.0 = IERR 

call IUNI (NMAX.NN»ALFT*NTAB» 
/AA3» I ORDER* alt, AR3. IPT3* I ERR) 
IERL(3,K )=I£RR 
DO 107 J=‘ ,8 


AZl ( J,K )=AR1 ( J ) 
AZ2(J*K)=AR2(J) 
107 AZ3( J*<)=AR3(J) 
105 CONTINUE 


IPT4=IPT5=-1 


CALL IUNI <N.VIAX,NN*ALr T ,NTAB, 

/Cl , I order. Alt* cc. ,ipta* ?e;?ci ) 
call IUNI (NyAX.NN.ALPT *NTAb* 
/C2. lORDER.Al T.CZ2* IPT5. IERC2 ) 
139 CONTINUE 


NMAX=NN=8 


NTAB=3 

determines the combination OF parameters relevant 
TO the case and sends control to appropriate 


SEGME 

NT 

Or 

the 

PROG 

,RAM 


if 

(OE 

• EG 

.0. 

) GO 

TO 

108 

IF 

(DW 

• EG 

• c. 

) GO 

TO 

109 

IF 

(OT 

• EQ 

.0. 

5 GO 

TO 

1 1 7 

GO 

TO 

1 1 1 





U IF 

(OW 

• EQ 

. I • 

) GO 

TO 

140 

IF 

(OT 

• EO 

.0 • 

> GO 

TO 

128 


GO TO 116 

140 IF (DT.e.0.0.) GO TO las 
K = 3 
OX = OT 
GO TO 130 
109 CONTINUE 

C COMES HE«E IF IT IS AN EMITTANCE VARIATION 
C CASE IN THE ORY ATMOSPHERE 
DO 112 N=1 « 3 
DO 112 J=1 »8 
AX( J.N)=AD( J»N) 

112 CONTINUE 

DO 129 N*1 .3 
129 A(N)=AX(3.N) 

IT=0 
DX = DE 

115 IT=IT+1 

IF ( IT.GT.IO) GO TO 113 
COR=A ( 1 )*DX+A <2 }*0> *DX 
ASTN=E3T-COR 
DIF=AbS l ASTN-ASTM ) 

IF (DIF.I.E.0,01 ) GO TO 114 

ASTM=ASTN 

IPT=- 1 

call IuN I (NMAX .Nn « TEM3 .NTAB . 

/AX. lORDER.ASTN. A . IPT, I ERR) 

1ERX( I T ) = lERR 
GO TO 115 

C COMES HERE IF IT IS A WATER vAPOR ONLY CASE 
128 CONTINUE 
K=2 

DO 131 J=1 .8 
AX( J. 1 )=AZI ( J.< ) 

AX ( J.2 ) = A22 1 J.K ) 

AX< J.3 )=A23 <J.< ) 

131 CONTINUE 

DO 132 N=1 .3 

132 A(N)=AX(3.N) 

I r ^o 

DX = D'-' 

133 lT=n 

IF ) GO TO 1 1 3 

COR = .. A { 2 ) »DX *nx 4-A ( 3 ) ♦OX » Dx »D\ 

ASTN = E:.T- '.OR 




oif=abs<astn-astm ) 

IF (DIF.LE.0.01 ) GO TO 114 

ASTM=ASTN 

1PT=-1 


call IUNI (NMAX»NN»TEMS»NTAB» 

/^AX« I0RDER»ASTN,A, IPT.IERR) 

1ERX< IT )=IERR 
GO TO 133 

COMES HERE IF IT IS AN EMITTANCE 
variation CASE IN THE WET ATMOSPHERE 
138 K=1 

DX“OE 

COMES here IF IT IS A TEMPERATURE PROFILE VARIATION 
CASE IN presence OF STANDARD WATER vAPQR BURDEN 
130 continue 

DO 134 J=l.e 
AX(J« 1 )=A2l <J»K ) 

AXlJ.a )=AZ2( J.< ) 

AX( J.3 )=AZ3( JtX ) 

AW(J« 1 )=AZ1 ( J.2 ) 

AW(J.2 )=AZ2< J»2 ) 

AW(J.3)=AZ3<J«2 ) 

Un CONTINUE 

DO 135 N=l .3 
A (N» =AX <3 .N ) 

135 C(N)=Ajtf(3.N) 

1T=0 

136 IT=IT+1 

IF ( I T.GT . 1 O ) GO TO 113 
CORX=A ( 1 )*OX+A (2 )*OX*DX 
CORW=C( 1 )+C(2)+C(3) 

COR=CORX+CORW 

astn=ebt-cor 

DIFrABS (ASTN-ASTm ) 

IF (DIF.LE.0.01) GO ~0 114 

ASTM=ASTN 

|PT1= IPT2=-1 

call IUNI (NMAX.NN. TEMp«NTAB. 

/Ax. lORDER. ASTN.A, IPTl , I ERR) 


IERX( IT )=!ERR 

CALL IUNI (NVAX.NN.TEMS.NTAB. 
/An/. ICRDER.ASTN.C, IPT2. lERR) 
IERW( I T )= lERR 
GO TO 1 36 



C COMES HERE IF IT IS A SURFACE EMITTANCE - 
C TEMPERATURE PROFILE VARIATION CASE IN THE 
C PRESENCE OF STANDARD WATER VAPOR BURDEN 
111 CONTINUE 

DO 118 Jst«8 
AX( Jt 1 )=AZ1 ( J« 1 > 

AX( J»E )3AZ2 ( J* 1 ) 

AX( J«3 )3AZ3( J« 1 > 

AW< J« 1 )=AZ1 (J«2 ) 

Aw( J-<2 ) -AZ2 ( J«2 ) 

AW( J«3 )=AZ3(J.2 ) 

AY( J. 1 )=AZ1 ( J.3 ) 

AY(Jt2 »®AZ2( Jt3) 

AY(J.3 )=AZ3( J«3 ) 

118 CONTINUE 

DO 119 N=1 .3 
A(N)=AX(3.N) 

8(N)*AY(3.N) 

C(N)=Aw(3.N) 

119 CONTINUE 
IT=0 

DX = DE 
OY=DT 

120 IT=n + l 

IF ( IT.GT.IO) GO TO 113 
CORX=A( 1 )*DX+A{2)*DX*DX 
CORY=3( 1 )-»-0Y 
CORW=C (1 )+C (2 )+C ( 3 ) 

cor=cokx+cory+corw 

AS7N=EbT-C0R 
D I F=A8S ( ASTN-ASTM ) 

IF (DIF.LE*O.0l ) GO TO 114 

ASTM=ASTN 

IPTl = IPT2= IPT3=-1 

CALL lUNI (NMAX*NN»TEMS»NTA o, 

/AX« IORDcW»ASTN«A. IPTl « I ERR) 

IERX( IT )=IERR 

call IUNI (NMAXf NN*Tc.vS.NTAH» 

/Ay» lOROER.ASTN.B, IPT2* I ERR ) 
iery( it ) = IERR 

CALL lUNI (NMAX.NNtTEMS.NTAB* 

XAw. lORDER.ASTN.C. IPT3. I ERR ) 

IERW( IT )= lERR 
GO TO 120 


(jHlOlisAL 
(■jir po<)F 


C COMES HERE !F IT IS A SURFACE EMITTANCE - 
C WATER VAPOR BURDEN VARIATION CASE BUT GOES TO 
C STATEMENT 138 IF THE WATER VAPOR BURDEN IS STANDARD 
117 CONTINUE 

IF (DW.EQ.l.) GO TO 138 
DO 137 N=l,3 
DO 137 J=l,8 
137 AX< J»N )=AD( J*N' 

DO 121 J-1 *8 
AY(J« 1 )=AZ1 ( J.2 ) 

AY( Jt2 )=A22 ( J.2 ) 

AY( J.3 )=A23( J.2 ) 

121 CONTINUE 
CK = C21 (3 ) 

DO 122 N=1.3 
A(N)=AX(3.N) 

B<N)=AY<3.N) 

122 CONTINUE 
1T=0 

DX = OE 
DY=DW 

123 IT=IT+1 

IF ( IT.GT.IO) GO TO 113 
CORX = A( I )*OX+A (2 ) *DX*DX 

CORY=B( 1 )*DY+B{2 ) #DY*OY+B ( 3 ) ^DY^Oy^O V 

C0R2 = C<*C0RX-»C0RY 

COR=CORXtCORY+CO»2 

ASTN=EoT-COR 

DIF = ABS ( ASTN-ASTM ) 

IF IDIF.LE.O.Ol ) GO TO llA 

ASTM=A3TN 

NTAB=3 

IPTl = 1PT2=IPT3=-1 

call lUNI (NMAX.NN.TEVS.NTAB. 

XAX. lORUER.ASTN.A, I^T I , I ERR ) 
lERX ( I T )= lERR 

call IUNI (NMAX.NN.TEMS.NTAB. 

/AY. IORDER.ASTN.B, IPT2. I ERR ) 
lERY ( IT )= lERR 
NTAB=1 

call I UN I (NMAX.NN.TEMS.NTAB. 

/C21 . 1 ORDER. ASTN.CK. IPT3. I ERR) 

IF.RC (IT) = IERR 
GO TO 1 23 


gUALrrY 


80 


u u 


COMES HER IF IT IS A WATER VAPOR BURDEN - 
TEMPERATURE PROFILE VARIATION CASE. 

I 16 CONTINUE 

DO 124 U=1 .8 
AX(J. 1 )-AZl ( Jf2 ) 

AX( J.2 )sAZ2 ( J«2 > 

AX ( J.3 ) sAZ3( J.2 ) 

AY(J. 1 )=AZl <U.3 ) 

AY( J.2 )sAZ2<U»3 ) 

AY( J.3 )=AZ3( J.3 ) 

124 CONTINUE 
C<=CZ2 (3 ) 

DO 125 N=l .3 
A(N)=AX(3«N) 

B(N)=AY<3»N) 

125 CONTINUE 
IT = 0 
DX*DW 
DY=DT 

126 IT=IT+I 

IF ( IT.GT.lO) GO TO 1 13 

CORX=A( 1 )*DX+A(2 )*DX*DX+A(3)*DX*DX*DX 
CORY=B( 1 )*OY 
CORW-A( 1 )+A(2)-t-A(3) 

CORZ = C<* C COR JK-CORX ) ♦CORY 
COR=CORX+CORY+CORZ 
ASTN=EBT-C0R 
DIF=ABS ( ASTN-ASTM ) 

IF (OIF.LE.0.01 ) GO TO 114 
ASTM=ASTN 
NTAB=3 

IPTI = IPT2= IPT3=-1 
CALL lUNI (NMAX.NN.TEMS.NTAB. 

/AX. IORDER.ASTN.A, IPTl , I ERR) 

IERX( IT ) = IERP 

CALL I UNI < NMAX .NN . TEMS .NTAB » 

/AY. lORDER. ASTN.B. IPT2 . I ERR ) 

IERY< IT)=IERR 
NTAB=1 

call I UN 1 (NMAX »NN « TEM3 .NTAB . 

/CZ2. IORDER. ASTN.C<» IPT3. lERP) 
lERC < I T ) = IcRR 
GO TO 126 

C RESETS THE CONVERGENCE PARAMETER Ir THE RESULT 


DID NOT CONVERGE AS 
LAST ESTIMATE EQUAL 

113 ICON=0 

114 surtemp=astn 

GO TO 127 
110 surtemp*ebt 
127 CONTINUE 

PRINT 60. SURTEMP 

62. DIF 

63. ICON 
64 
61 . 

61 . 

61 . 

61 • 

61 . 

61 . 


DESIRED. AND 
TO SURTEMP 


SETS the 


ORIGINAL FAGL IS 
OF POOR QUALITY 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


>3) 


60 

62 

63 
61 

64 


( ( lERLCN.X >.N=1 .3) .K=l 
lERCl .IERC2 
( lERX < I ) . I =1 . 10 ) 

( lERY ( I ) . I =1 . 10 ) 

( lERC (I ) . I =1 . 10 ) 

( lERW (I ) . I =1 . 10 ) 

FORMAT! IHl///. 15X.* SURFACE TEMPERATURE =*.Fa.2) 
FORMAT (//.15X.* I TEWA I 1 UN Hb^IUUAL 
FORMAT!//. 15X.* CONVERGENCE PARAMETER =*.17) 
format !//. 1 017 ) 

format !///. 12X.* ERROR PARAMETERS FOR INTERPOLATION*) 
STOP 


